Seebeck effect in dilute two-dimensional electron systems: temperature dependencies 
of diffusion and phonon-drag thermoelectric powers 
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Considering screeening of electron scattering interactions in terms of tlie finite-temperature STLS 
theory and solving the linearized Boltzmann equation (with no appeal to a relaxation time approxi- 
mation) , we present a theoretical analysis of the low-temperature Seebeck effect in two-dimensional 
semiconductors with dilute electron densities. We find that the temperature (T) dependencies of 
the diffusion and phonon-drag thermoelectric powers {Sd and Sg) can no longer be described by the 
conventional simple power-laws. As temperature increases, \Sd\/T decreases when T > O.lep (ef is 
the Fermi energy), while \Sg \ first increases and then falls, resulting a peak located at a temperature 
between Bloch-Griineisen temperature and ef- 
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I. INTRODUCTION 

Heat generation with increasing density of integrated 
electronic devices is one of the serious restrictions block- 
ing the further development of conventional electron- 
ics. To overcome this obstacle, there have been pro- 
posals to use heat to storage and transport information. 
Such an intriguing concept provides hope of discover- 
ing new physics and new green technology, stimulating a 
great deal of theoretical and experimental investigation. 
Recently, two new subfields associated with heat elec- 
tronics have emerged, namely "Phononics"^ and "Spin 
Caloritronics" Phononics is devoted to the use of heat 
current to perform computational operations and many 
heat devices, such as heat diodes^ heat transistors^, 
thermal memorj^ etc., have been proposed and/or con- 
structed. Spin Caloritronics concerns on the motion of 
magnetization or of electron spin induced by heat or a 
temperature gradient j2, In this field, many phenomena, 
such as a spin Seebeck effect thermodynamic control of 
magnetization in ferromagneto-nonmagnetic structuresji 
etc. have been reported. 

To use heat in electronics, the conventional method 
is usually based on thermoelectric effects, which con- 
vert a temperature gradient to electric voltage. Among 
them, the Seebeck effect (SE), initially discovered in 
metals by T. Seebeck in the 1820s, ^'^ is widely used 
for thermoelectric generation and for temperature sens- 
ing. The first observations of SE in bulk semiconduc- 
tors, such as Germanium^'^'^^ and Silicon,^? were re- 
ported in the 1950's. With the recent development 
of technology in the fabrication of semiconductor mi- 
crostructures, investigations of thermoelectric effects in 
two-dimensional electron gases (2DEG) have been car- 
ried out both experimentalljJ^"— and theoretically<^"— 

It is well known that, in thermoelectric power (TEP), 
which is the main characteristic quantity of SE, there are 
two components; namely, diffusion thermopower, Sd, and 
phonon-drag thermopower, Sg. At relatively low tem- 



perature, the diffusion process in the Seebeck effect has 
been expected to be dominant since the electron-phonon 
scattering is relatively weakJ^"— However, a careful 
analysis of experimental data indicates that phonon-drag 
in 2DEGs also plays an important role even at temper- 
ature T < 10K^iai20i36-4i^^^ Furthermore, there 
have also been studies of a sign change of diffusion TEP 
in a Si-MOSFET,i^ and of the effects of weak local- 
ization on TEP,^^-^^-i^ as well as the TEP of composite- 
fermions,2i'23.24.33,49^-^ ^^^^ oscillation of TEP in low 

magnetic field )^'''^^i^° etc. 

To understand the microscopic mechanisms in SE, it is 
necessary to separate Sd and Sg from the total TEP that 
is measured. One way to do this is to analyze the tem- 
perature dependencies of Sd and Sg. In the absence of 
phonon-phonon scattering in the phonon relaxation pro- 
cess, low-temperature Sg vs T behavior has been taken 
to be of the form: Sg oc AT" with A as the phonon mean 
free path and rt = 3 or n = 4 for dirty or clean sam- 
ples respectivelyi^Si^ The diffusion TEP has often been 
assumed to vary linearly with temperature.™ However, 
Sankcshwar, et al. showed that the inelastic feature of 
electron-phonon scattering may result in a nonlinear tem- 
perature dependence of Sd in relatively clean 2D samples 
in the Bloch-Griineisen (BG) regime, i.e. T < Tbg [the 
BG temperature Tbg = '^kpUsX with Ugx as the phonon 
velocity in branch A is about 5 K for a 2DEG with typical 
density ~ 10"'^"'^ cm~^]i^ A few experiments were de- 
voted to the direct measurement of diffusion TEP .^^'^^i'^^ 
Ying, et al. observed pure diffusion TEP only in the case 
T < 0.2Ki2i Recently, using hot electron thermocouple 
structures, the diffusion TEP has been directly detected 
by Chickering, et al. when T < 2 

It should be noted that the simple power-laws of Sd 
and Sg vs T, obtained previously, were derived on the 
basis of a relaxation time approximation (RTA), which 
is valid when T, Tbg ^ Ef- Recently, motivated by the 
observation of a so-called metal-insulator transition in re- 
sistivity vs temperature, clean undoped heterojunctions 
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with electron density Ug as low as Ug < 10"'^'' cni~^ have 
been studied extensivelyi^""— In these systems, T and 
Tbg are comparable with the Fermi energy even at low 
temperature and therefore deviations of Sd and Sg vs 
T from the conventional results are expected to be ob- 
served. 

In this paper, within the framework of Boltzmann 
equation, we present a theoretical investigation on ther- 
moelectric effects in 2D electron GaAs/AlGaAs systems 
with carrier densities rig = 0.23 ~ 1.06 x 10^°cm~^. To 
account for the screening of scattering interactions in a 
2DEG with such low Us, the finite-temperature Singwi- 
Tosi-Land-Sjolander (STLS) theory, a scheme beyond 
random phase approximation (RPA), is employed.— 
Furthermore, to carefully treat inelastic electron-phonon 
scattering, the Boltzmann equation is solved with no ap- 
peal to a relaxation time approximation, using an en- 
ergy expansion method. Das Sarma and Hwang have al- 
ready presented a qualitative explanation of experimen- 
tal observations of resistivity in a 2DEG with such dilute 
rig by means of a Boltzmann equation combined with 
RPA-screened electron-impurity scattering. ^''"^^ In the 
present paper, performing numerical calculations with 
STLS screening appropriate to low carrier density, we 
find that the temperature dependencies of Sd and Sg in 
dilute 2D systems are significantly different from those 
in the high-electron-density limit. When temperature 
increases, \Sd\/T no longer remains unchanged: it de- 
creases for T > QAep- In our calculation of phonon-drag 
TEP vs temperature, a peak appears: \Sg\ first increases 
and then falls as temperature increases. 

The paper is organized as follows. In Sec. II, an en- 
ergy expansion method for solving the Boltzmann equa- 
tion beyond the RTA is presented along with the self- 
consistent finite-temperature STLS theory. Numerical 
investigation of the temperature dependencies of diffu- 
sion and phonon-drag TEPs for various dilute electron 
densities are exhibited in Sec. III. Our results and con- 
clusions are summarized in Sec. IV. In the Appendix, we 
also provide analytical results for Sd and Sg vs T in the 
high-electron-density limit, obtained by the energy ex- 
pansion method. 



II. THEORETICAL CONSIDERATIONS 
A. Electron and phonon Boltzmann equations 

When a two-dimensional electron with momentum p = 

2 

(pcos(/?p,psin</3p) and energy Sp — (to* is the effec- 
tive electron mass) is subjected to a weak electric field 
E and a thermal gradient VT, its kinetic motion can be 
described in terms of a nonequilibrium distribution func- 
tion, /p, which is determined by a linearized Boltzmann 
equation of form 



eE ■ 



T 



dfojsp) 



(1) 



Here, /z is the chemical potential, T is the lattice temper- 
ature, fo{ep) = {exp[(ep — //)/T]-|-l}~^ is the equilibrium 

p/to* 



electron distribution function and Vr 



VpEp 



is the electron velocity. In Eq. ([T|) , /scatt is the scatter- 
ing term due to electron-impurity and electron-phonon 
interactions and it can be written as /scatt — /imp + /ph. 
/imp represents the contribution to /scatt from electron- 
impurity scattering: 



/imp = ^27r}^\Ug\'S{ep - £p_q) (/p - /p_q) , (2) 



while /ph is associated with the electron-phonon interac- 
tion: 

/ph = -27r \MQxf6{ep±q-epT^Qx) 

Q,A,± 

X [^qa/p(1 - fp±i) - N^xfp±ii^ - fp)] ■ (3) 

In Eqs. (HI) and ([3]), Uq is the electron- impurity scat- 
tering potential, Mqx is the matrix element for inter- 
action between the 2D electrons and 3D phonons, and 
N^^ = Nqx + ^ T 5- ^QX and Nqx, respectively, are 
the energy and number of nonequilibrium phonons with 
three-dimensional momentum Q = {q,qz) = i<lx,Qy,<lz) 
in branch A. 

Since the temperature gradient may drive the phonons 
out of equilibrium, Nqx in Eq. ([3]) differs from the num- 
ber of equilibrium phonons, uqx = [exp{Q,Qx/T) — 1] , 
and it is determined by the Boltzmann equation for 
phonons: 
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Here, ^ ^^^^ ^ is the drift term, taking the form 



fdN 



QX 



V dt 



-Uqx ■ VNqx, 



(5) 



with Uqx as the phonon velocity. Note that, in the 
present paper, the magnitudes of uqx are assumed to 
be independent of Q and they are denoted by Usx (longi- 
tudinal and transverse acoustic phonons are denoted by 

Usi and Ust, respectively). {—^^^ is the relaxation 

term due to the boundary and phonon-phonon scatter- 
ings, written as 
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{Nqx ~ Uqx) 



(6) 



with Tbs as the relaxation time due to boundary scatter- 
ing, l/rbs — Uqx/ h-r^ and Tpp as the relaxation time due 



to phonon-phonon scattering, l/rp 



AxT^^x^ 
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{~W^] '^^ phonon scattering rate due to the 

V / ep 



3 



electron-phonon interaction, as given by 



dNo\\ 27r , — , „ ^ 

^ ^ =-— 5s|MQA|2^{±<5(£p~£p±q±r!QA) 
cp p.± 



dt 



X TV, 



(7) 



with gs as the spin degeneracy and L as the sample size 
along the direction perpendicular to the 2D sheet. In 
the case of a weak temperature gradient, Eq. ([5]) can 
be solved analytically and the steady-state number of 
nonequilibrium phonons can be written as 



nqx ~ TpxUQx ■ VT 



dfii 



QA 



dT 



(8) 



Here, 1 / Tp\ = 1/rbs + 1/ Tpp + 1/ Top and Top is the phonon 
relaxation time due to electron-phonon scattering, taking 
the form 
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^ E \MQxfS{ep-g-ep+nQx)[foiep-g)-foisp)] 
p 

(9) 



Energy-expansion method to solve electron 
Boltzmann equation 



To solve the electron Boltzmann equation, Eq. ([1]), we 
assume that the nonequilibrium distribution function fp 
takes the form 



fp = fa{£p) +9p 



dfaisp) 



de 



(10) 



with gp as an unknown function. In previous studies, 
when electron-optical-phonon scattering can be ignored 
at low temperature, gp is usually obtained using the 
relaxation time approximation (RTA). Obviously, RTA 
is valid only in the high-electron-density limit. In the 
present paper, in order to study diffusion and phonon- 
drag TEPs for relatively low electron density, we follow 
the idea proposed by Allen for an investigation of trans- 
port in metalsj^ which assumes that gp can be expanded 
in terms of basis functions XJn{p)' 



9p 
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CjnXJn (P), 



(11) 



with Cjn as the coefficients of expansion. In a 2D system 
with a parabolic dispersion relation, the functions XJn{p) 
can be written as 



XJn{p) = 47rF,7(p) 



■nn{£p) 



(12) 



Here, Fj{p) are the basis functions for the expansion of 
gp with respect to the angle of the momentum vector p, 
and they can be chosen as sine or cosine functions of mul- 
tiples of the angle (fp. rin{s) are n-th order polynomials 



in electron energy e and they are orthogonal with respect 
to the weight function — 9/o(e)/9e: 

(-^^) Vnie)vrn{e)de = <5„,„. (13) 

It is noted that to study the transport in metals, the lower 
limit of energy integration in Eq. (1131) can be assumed to 
be —oo, since the Fermi energy in metals usually is much 
larger than the bottom of electron energy bandj^ How- 
ever, in three- or two-dimensional semiconductors, the 
finite bottom of the energy band or subband may affect 
transport properties, especially at relatively high tem- 
perature (or for relatively low electron density). Hence, 
in Eq. (1131) . the lower limit of integration is maintained 
equal to zero. Further, in our study, we assume that 
77„(e) take the form 



Vnis) 



m=0 



(14) 



which also differs from that proposed by Allen. In 
Eq. (jl4p . the parameters anm are determined from the 
orthonormality conditions of r]„{e). In general, they are 
independent of s but may depend on the lattice tem- 
perature, as well as on the Fermi energy fi. Note that 
for n — 0, 770(e) is an energy- independent constant: 
Vois) = Vo- 

Furthermore, without loss of generality, we assume 
that the electric field and temperature gradient are ap- 
plied along the x axis. Thus, in 2D semiconductors with 
parabolic dispersion, only one term with basis function 
Fj^x (p) — cos fp need be considered in the expansion of 
gp with respect to ipp. Multiplying both sides of Eq. ^ 
by Xn{p) [Xn{p) = Xn,j=x{p) = 'iiT COS (ppT]n{ep) / p] and 
performing the summation over p, the linearized Boltz- 
mann equation for electrons can be rewritten as 



eE 



T 



m—O 



i7m H —D„ 



T 



n'=0 



Qnn' Cji' 



with 7„ 



-m 9/o(ep) 



(15) 

In this equation, the 



Jo de{e-fi)e' —^^ 
third term on left-hand side is the source of the phonon- 
drag effect; it describes the interaction between equilib- 
rium electrons and nonequilibrium phonons. In it, £)„ 
take the form 

Dr^ = 27r ^ (=Fl)|MQApx„(ep)(5(ep_q -SpT ^QxVpX 

Q,P,X,± 



xu. 



^QX 

QX rp 



"-qa(1 + nqx) [fo{£p) - /o(£p-q)] , (16) 



with Uq^ as the x component of uqx- Note that Eq. 
is derived from Eq. ([3]) by substituting into it the ex- 
plicit form of the number of nonequilibrium phonons, i.e. 
Eq. ([5]). On right-hand side (r.h.s.) of Eq. Qnn' are 
associated with the scattering term Isc and they can be 
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written as Q„„/ = Q™? + Q^^, with Q™? and QP^,, re- 
spectively, taking the forms {fpq is the angle between p 



and q) 



dfoisp) 



X(5(£p - £p-q) COS (yCpq 



(17) 



and 



J2 \MQ,x\''Siep-q-epT^Q,x> 



p,Q,\,± 



x/o(ep)[l - /o(ep-q)] — 



{p-qcosippq) 



Zp-qV^r 



(18) 



Thus, the original linearized Boltzmann equation is re- 
duced to Eq. (fT5|) . a system of linear equations for C„. 
After C„ are determined, the macroscopic charge current 
can be evaluated through 



p 



dfoisp 



9pVp 



-gscCo/vl (19) 



Since there are three driving terms in Eq. (jisp . its so- 
lution, Co, can be written as Co = Cq 



(c) 



C, 



(d) 



c, 



(g) 



with C^'^\ C^'^\ and C^^^ determined from Eq. (HH) in 
the presence of only the first, the second, or the third 

(c) 

driving term, respectively. Obviously, Cq is propor- 
tional to E and it determines the conductivity as cr = 
-gse&Q^ /{-qlE). and are proportional to V^T 
and they are associated with the diffusion and phonon- 
drag TEPs, respectively: Sd — —g8eC'^^/{rjl'S/xTa) and 

We note that such an energy expansion method pre- 
sented here can also reproduce the previous RTA results 
in high-electron-density limit. We present a detailed cal- 
culation of the high-ris TEP as a function of T in the 
Appendix, considering screened electron-impurity scat- 
tering as well as screened piezoelectric interaction and 
unscreeened deformation interaction between electron- 
acoustic phonons. There, the well-known Mott relation 
is obtained for Sd^ and the lowest-order correction to the 
Mott formula at low temperature may come not only 
from electron-impurity scattering but also from the in- 
teraction between electrons and acoustic phonons when 
T lies within the equipartition (EP) regime, Tbg -^T <^ 
ep- We also obtain the well-known law for Sg vs 
temperature for T within the BG regime. Besides, in 
the presence of only boundary scattering in the phonon 
relaxation process, Sg is found to be independent of tem- 
perature when Tbg ^ T <C Ef- 



C. Finite-temperature STLS theory 

To analyze resistivity as a function of T in dilute 2D 
systems, it is necessary to clarify the role of screening in 
electron-impurity and electron-acoustic-phonon scatter- 
ings. Using RPA-screened electron-impurity scattering. 
Das Sarma and Hwang have qualitatively explained the 
experimental observations in a dilute 2DEG.— How- 
ever, in the GaAs systems that we study, the dimen- 
sionless Wigner-Seitz density (or interaction) parameter 
Ts = l/(aB-y/7m7) [aB = Atteok/ {m*e^) is the effective 
semiconductor Bohr radius and k is the background di- 
electric constant] can reach the value ~ 11.6 for 2D GaAs 
with Us — 0.23 X 10~^°cm~^. Hence, the local-field cor- 
rection to RPA is quite important and we therefore use 
the finite-temperature self-consistent STLS theory here. 

Within the framework of the Boltzmann equation ap- 
proach with interaction screening included, the scatter- 
ing potential is usually divided by the dielectric function 
e{q,uj), which takes the form 

e{q,L,) = 1 - VgH{q)[l - Giq)]xo{q,oj). 

Here, xoilT'-^) is the density-density correlation func- 
tion of the free 2D system, H{q) is the form factor of 
the electron-electron interaction in the 2D system, and 
Vq — e^/{2neQq) is the 2D Coulomb potential. G{q) is 
the static local-field factor whose value depends on the 
approximation that used. In RPA, G{q) is zero, while 

G{q) = q/ ^2-^/(72 + in Hubbard's approximationi2£ 

In STLS theory which we use here, the local field factor 
is determined by the structure factor S{q) through 



Giq) 



(20) 



On the other hand, S{q) is also related to G{q) via 
T 



Siq) 



E x(9,2i7rnT), 



(21) 



n— — 00 



with x(q, iw) = Xoilii^) / ^ilT^^) ^-s the response func- 
tion. Thus, Eqs. (|20|) and ([21]) form a closed system of 
equations, to be solved self-consistently by iteration. 



III. RESULTS AND DISCUSSION 

We carry out numerical calculations to investigate the 
thermoelectric effect of a dilute 2D electron gas in a 
GaAs/AlGaAs heterojunction at temperature T < 5K. 
The electron Boltzmann equation is solved by means of 
the energy expansion method and the screening of scat- 
tering is evaluated self-consistently within the framework 
of the finite-temperature STLS theory. In these calcula- 
tions, the screened electron-impurity scatterings due to 
both remote and background impurities are considered. 
The corresponding scattering potential takes the formal 



PqWHqM 



(22) 
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with = V^lN.e-^i'Iiqf + N^J{q)/q], I{q) and 

J{q) are the form factors, N\, is the density of background 
impurities, and represents the density of remote im- 
purities focated at distance s from the heterojunction 
interface on the AlGaAs side. 

In regard to the electron-phonon interaction, only 
acoustic phonons contribute to scattering at low tem- 
perature. The corresponding potential can be written 



IM, 



QM 



IM. 



QA| 



(23) 



with |MqaP as the matrix element of the electron- 
phonon interaction in three-dimensional plane- wave rep- 
resentation. In present paper, we consider both the defor- 
mation and piezoelectric interactions between electrons 
and acoustic phonons. It is well known that only the 
longitudinal acoustic phonon (LA) mode gives rise to de- 
formation scattering with matrix element 



lAf, 



QiAldcf 



2du^ 



(24) 



Here, d is the mass density of crystal and S is the shift 
of the band edge per unit dilation. Both the longitudinal 
and transverse (TA) acoustic phonons contribute to the 
piezoelectric interaction. The corresponding scattering 
matrix elements take the forms^ 



lAf, 



2 _ 3277^6^6^4 9ggg^g^ 



and 



11^ |2 ^ 327r^6^6i4 



with ei4 as the piezoelectric constant. 

In Eqs. ([24|) and ([25|) . the unscreened form of the 
electron-acoustic-phonon scattering through deformation 
potential is used, while the piezoelectric interaction is as- 
sumed to be dynamically screened. Such a treatment is 
based on the fact that these two interactions have com- 
pletely different origins. It is well known that piezoelec- 
tric electron-phonon scattering comes from the Coulomb 
interaction of electrons in an electric field induced by 
thermal vibration of atoms, and hence it is effectively 
screened by electron-electron interactions. However, the 
deformation scattering mainly results from the overlap of 
electron wave functions between different atoms in dis- 
torted lattices.— Thus, the deformation interaction be- 
tween electrons and phonons does not directly relate to 
the Coulomb interaction, and therefore it is inappropri- 
ate to use the screened form for the deformation poten- 
tial. We note that, employing the unscreened form for 
the deformation potential with the appropriate parame- 
ter S, good agreement between theory and experiments 



has been reached in a previous study on phonon-drag 
thermoelectric effect.'*^ 

In our numerical calculations, the parameters are cho- 
sen as follows: k — 12.9, d ~ 5.31 g/cm'^, Usi — 5.29 x 
10^ m/s, Ust = 2.48 x 10^ m/s, S = 8.5 eV, m* = 0.067too 
(too is free electron mass), ei4 — 1.41 x 10^ V/m. Since 
we are interested in the temperature and electron-density 
dependencies of the diffusion and phonon-drag TEPs at 
low temperature (T < 5K), the relaxation of phonons 
due to phonon-phonon scattering can be ignored and 
only the temperature-independent boundary scattering 
need be considered. Furthermore, the phonon mean free 
path is assumed to be A = 2.42 mmi^ The truncation 
of summation in the expansion of g„ is estimated by the 
convergence of the numerical scheme. We find that, for 
Us > 0.2 x 10^"cm~^ and T < 5K, rimax = 4 is sufficient 
to reach the required numerical accuracy. 

The low-temperature transport properties depend sen- 
sitively on the impurity densities. In the present paper, 
to obtain results in qualitative agreement with the exper- 
imental resistivity data of Ref. GOj, the density of charges 
in the depletion layer is chosen to be A^dep = 7 x 10^ cm~^ 
and the background impurity density is assumed to be 
constant: Nh = 1 x 10^^ m~'^. The remote impurity den- 
sity, N^, is determined from the mobility at T = 30 mK 
by assuming s = 210 nm. In Fig. 1, we plot the depen- 
dencies of resistivity p = l/cr on temperature for various 
electron densities. An evident "metal-insulator" transi- 
tion can be observed: when T increases p increases for 
dense n^, while it decreases for dilute rig. Such behavior 
of p vs r almost agrees quantitatively with experimental 
data in the case 1 K < T < 5 K for all Us which were 
studied (see Fig. 2 in Ref. 1 60). However, in Fig. 1, we do 
not see the small peaks for intermediate Us, which have 
been observed experimentally.— This is associated with 
the fact that the observed small peaks in p vs T are the 
result of weak (or strong) localization, which is ignored 
in our study. 

In Fig. 2, we plot the temperature dependencies of 1 5^1 
and \Sd\/T in a 2D GaAs heterojunction for various elec- 
tron densities in the range 0.23 < < 1.06 x lO^^cm"^. 
From Fig. 2(a), we see that, with an increase of tem- 
perature, \Sd\ increases. However, this increase is no 
longer linear. To clearly show the nonlinear dependen- 
cies of Sd on T, the temperature dependencies of \Sd\/T 
for r < 5K are plotted in Fig. 2(b). We see that when 
temperature increases, \Sd\/T remains constant only for 
T < O.lep- Beyond this regime, \Sd\/T decreases with 
an increase of temperature. Such nonlinear dependence 
of Sd on temperature mainly stems from broadening of 
the electron distribution function at relatively high tem- 
perature. 

It should be noted that electron-phonon scattering also 
may affect \Sd\ at relatively high temperature. To show 
this, in Fig. 3, S'd vs T is plotted both in the absence and 
in the presence of electron-phonon interactions. It is clear 
that the contribution from electron-phonon scattering to 
Sd is important for relatively high n^. This is associated 
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FIG. 1. (Color online) Temperature dependencies of resistiv- 
ity for various densities of electrons (from the top): Us — 0.23, 
0.29, 0.36, 0.42, 0.55, 0.68, 0.80, and 1.06 x 10^° cm~^ The 
remote impurity sheet is assumed to be located at 210 nm 
from the interface on AlGaAs side. The density of charges in 
the depletion layer is assumed to be A'dep = 7 x lO'' cm~^ 
and a constant density of background impurities is used: 
iVb = 1 X 10^* m-^ 



with the fact that for dilute n^, electron-impurity scat- 
tering is so strong that the electron-phonon interaction 
is relatively unimportant within the temperature regime 
studied. From Fig. 3, we also see that the magnitude of 
Sd in the presence of electron-phonon interaction is al- 
ways less than that in the absence of electron-phonon 
scattering, reflecting the fact that contribution to \Sd\ 
from electron-phonon scattering is negative. 

In Fig. 4, we plot the temperature dependencies of 
phonon-drag thermoelectric power for various electron 
densities. We see that Sg vs T for dilute Ug is signifi- 
cantly different from that in 2D systems having a dense 
electron density. It is clear that for relatively high rig 



(for example in the case Us 



1.06 X lOiOcm-2) 



creases as T increases and then it saturates at a relatively 
high temperature. This can be explained qualitatively by 
means of the asymptotic behavior of Sg in high-rij limit, 
presented in the Appendix: when T 'ti sp, Sg first in- 
creases with an increase of temperature as ~ and it 
becomes independent of temperature at high tempera- 
ture. From Fig. 4, we also see that for relatively low n^, 
a peak appears in 15*31 vs T. The position of peak de- 
pends on electron density: the peak moves towards the 
low temperature side with a decrease of electron density, 
but it lies always between Tbg and ep- 

It should be noted that the appearance of a peak in Sg 
vs T can be understood as the result of competition be- 
tween (i) broadening of the Fermi distribution function 
and (ii) decrease of the rate of nonequilibrium phonon 
production, induced by an increase of temperature. As T 
increases, the Fermi distribution broadens and dragging 




r(K) 

FIG. 2. (Color online) Temperature dependencies of (a) \Sd\ 
and (b) \Sd\/T for various densities of electrons (from the 
top): Us = 0.23, 0.29, 0.36, 0.42, 0.55, 0.68, 0.80, and 1.06 x 
lO^^cm"^. The dashed line in Fig. 2(a) indicates the Fermi 
energy ef (in Kelvin). 



electrons out of equilibrium by nonequilibrium phonons 
is facilitated. As a result, the phonon-drag TEP \Sg\ 
increases with increasing T. However, as T further in- 
creases, the rate of nonequilibrium phonon production 
induced by a temperature gradient decreases, leading to 
a decrease of Sg with the further increase of T. Com- 
petition of these two factors results in the nonmonotonic 
dependence of Sg on T. 

From analysis presented above, it is clear that to ob- 
serve the nonmonotonic dependence of jSgl on T in the 
presence of only boundary scattering in the phonon relax- 
ation process, two conditions are required. One condition 
is that the Fermi energy should be much smaller than the 
critical temperature at which phonon-phonon scattering 
is important in phonon-relaxation. In 2D GaAs systems, 
such a critical temperature is about 10 K,^*' leading to an 
estimate of electron density in a 2D GaAs system for ob- 
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r(K) 

FIG. 3. (Color online) Effect of electron-phonon scattering 
on diffusion TEP for Us = 0.23, 0.42, and 1.06 x lO^cm"^. 
The solid and dotted lines indicate \Sd\ in the presence and 
in the absence of electron-phonon scattering, respectively. 

servation of the peak in Sg vs T as < 2.5 x 10^" cm^^. 
The second condition for observation of the peak in Sg vs 
T is that Tbg should be comparable with sp- Otherwise, 
the peak disappears: as temperature increases, Sg mono- 
tonically increases and reaches saturation at a relatively 
high temperature when Tbg <SC £_f. 

Note that in typical thermoelectric experiments, the 
measurable quantity is the total thermoelectric power, 
causing difficulty to separate the diffusion and phonon- 
drag contributions. In the dilute 2D systems studied 
here, the phonon-drag TEP is dominant over almost the 
whole temperature regime T < 5 K. For example, for 
Hs = 1.06X 10i"cm-2, \Sg\ exceeds \Sd\ when T > O.IK. 
Hence, in the temperature dependence of the total TEP 
S = Sd + Sg, which is plotted in Fig. 5, the features are 
almost the same as those in Sg vs T. Therefore, to ob- 
serve the nonlinear dependence of diffusion TEP on T as 
shown in Fig. 2, specific structures, such as a hot-electron 
thermocouple;^ are required. 



IV. CONCLUSIONS 

Employing the energy expansion method to solve the 
Boltzmann equation and taking account of the screen- 
ing of interactions in terms of finite-temperature STLS 
theory, we have carried out a theoretical investigation 
of the thermoelectric effect in a two-dimensional elec- 
tron system with dilute electron density 0.23 < n-g < 
1.06 X lO^^cm"^. The temperature dependencies of the 
diffusion and phonon-drag thermoelectric powers have 
been carefully analyzed for T < 5 K and our results ex- 
hibit deviations from the conventional simple power laws. 
We find that, in dilute 2D systems, \Sd\/T remains con- 



T(K) 

1 2 3 4 5 




r(K) 

FIG. 4. (Color online) Temperature dependencies of (a) \Sg\ 
and (b) \Sg\/T* for various densities of electrons (from the 
top); Us = 0.23, 0.29, 0.36, 0.42, 0.55, 0.68, 0.80, and 1.06 x 
lO^'^cm"^. The dotted and dashed lines in Fig. 4(a) indicate 
Tbg and ef (in Kelvin), respectively. 



stant only when T < O.lep and it decreases with an 
increase of temperature out of this regime. We also 
observe a peak in the temperature dependence of \Sg\, 
which arises from competition between thermal broad- 
ening of distribution functions and decrease of the rate 
of nonequilibrium phonon production, induced by a tem- 
perature increase. 
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r(K) 



FIG. 5. (Color online) Temperature dependencies of magni- 
tude of total thermoelectric power S = Sd + Sg for various 
densities of electrons (from the top): Us — 0.23, 0.29, 0.36, 
0.42, 0.55, 0.68, 0.80, and 1.06 x 10^° cm'^ 



Appendix: Asymptotic behaviors of 5*^ and 5*3 in the 
high-electron-density limit 



We verify that the energy-expansion method presented 
in section II produces the conventional expressions for 5*^ 
and Sg in the case T ^ ep [sf = STrn^/ (^s'ti*) is the 
Fermi energy]. Obviously, in high-electron-density limit, 
it is sufficient to consider only the two lowest terms in 
the expansion of 5p: gp ~ [Co?7o + Ci 771(e)] cos (/3p. Based 
on the orthonormality conditions of 77„(e), Eq. p^ . 770(e) 
and 771(e) can be written, respectively, as 



1 



Vim 



= 1 



-e 
2 



3 

-e 



-HO(e-^) (A.l) 



and 



7?i(e) 



V3{eF — e) J 3\/3(e_F — e) 
^ ^ I 27r3T 



1 /ep 



2 \ T 

-0 (t-^ 



^ + 1 
T 




(A.2) 



From Eqs. (|A.ip and (|A.2p it is obvious that the cor- 
rections to the leading terms of 770,1 (e) are exponen- 
tially small for T <C ep and hence r?o(e) = 1 and 
771 (e) = •\/3(£_F — ^)/{'^T) can be used in the calculation 
that follows. 



1. Temperature dependence of resistivity in the 
high-electron-density limit 

Before proceeding to analyze Sd and Sg vs T, it is 
useful to evaluate the temperature dependence of resis- 
tivity, which is defined as p — 1/cr = m* {t^^') /{use^) 
((r-i) is the average inverse relaxation time) and is 
proportional to Qoo in high-electron-density limit: p = 

4Qm/{9se') = TiliQZ" + Qt,)/{9se'). 

We first consider the temperature dependence of Qq™^ 
resulting from RPA-screened electron-impurity interac- 
tion [G{q) = is used]. Using the potential given by 
Eq. (22), Qoo'' can be expressed as 



^imp 
'00 



de5{eF - e)A™P(e), 



with A™P(eF) taking the form 
A™P(ef ) = ^ / dq 



£(9,0) 



(A.3) 



.(A.4) 



In the low-temperature limit, A""P(eF) can be further 
expanded as 



A'-P(eF) = ^ 



X(7t + 1) 



with (j){x) — dy 




Vv 

cosh*^ {y-\-x) 



eT=o(2fcF,0)]"+2 

/■°° dx 

/ ^0(x), (A.5) 
Jo 

12 From Eqs. (IA.3|) and 



(jA.sp it follows that, for T ^ e^, Qqo^ can be written as 
Q[,7 = A-P(eF)=g^ 



•)imp 
^00 



T=0 



with r™P(eF) determined by 

8777* V2kFH2kf 



'(sf) 



Tre% [eT=o{'2kF,0W Jo y/x 



r™p(eF)r + o(r3/2)^ 

(A.6) 



°° dx 

^<t>{x), (A.7) 



and Q'™^ 



'00 



T=0 



is obtained from Eq. (jA.5|) by setting T 



0. From Eq. (|A.6|) it is clear that the first-order finite- 
temperature correction to is linear in T, consistent 
with previous transport studies. ^^'^^ Note that such a 
correction comes mainly from the temperature depen- 
dence of the dielectric function in the screened electron- 
impurity scattering potential. Moreover, if the screening 
of electron-impurity scattering is considered by means of 
the finite-temperature STLS theory, an additional tem- 
perature dependence associated with G{q) needs to be 
taken into account. 

Further, considering both the deformation and piezo- 
electric scatterings, we carry out the determination of the 
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temperature dependence of Qqq both in the BG (T ^ 
Tbg ^ £f) and m the equipartition (Tbg 'ti T <ti ep) 
regimes. In both cases, Tbg ^ £f and hence we can 
make the approximation: 

/o(ep)[l - /o(ep + ^Qx)] ~ (1 + nQx)^QxS{ep - ep)- 

(A.8) 

Performing the p- integration in Eq. (jl8p , Q^^/ for ?i = 
n' = can be written as 

QP^ ^ y de(S(ei. - £)AP'^(£) (A.9) 
with AP'^(e) defined by {kp = \/2m*e^) 

x^^QAG±(g, l^QA)nQA(rtQA + 1). (A.fO) 

Here, G±{q,nQx) = 2(£q ± nQx)[{kFqlm*f - {^qx ± 
eq)^]~^/^ is associated with (/Spg-integration over 6 func- 
tion. 

In the BG regime, T < Tbg < £f, G±{q,nQx) « 
2TO*(e, ± f^QA)/(fci.9), 1, and e((j,rjQA) « 

Qs/g with = m*e'^ /{2iTeoK) as the screening wave 
vector. Thus, Eq. (|A.10p can be rewritten in the low- 
temperature hmit as 



rp4 



dqz 



X / dq\MQx\ 
'0 



X J -oo 



1 



(A.ll) 



with Q — ^Jq^^\^(^ . Substituting the exphcit form of 
Mqx into Eq. ()A.11| . the contribution to Qqq in the BG 
regime from the deformation potential, Qcxf is given 
by 

T^S^ 

'00 " ^ 



Sndkpepi 
15C(5)T5s2 



2dkF£pU 



2 „,6 ' 



1)' 



(A.12) 



and the contribution from the screened piezoelectric in- 

, can be written as [Bla = 



teraction to < 



^ph-BG ^PZ-BG 

'Oa 



5 00 

— (7^ 



tjqf' and Bta = {(i 
^PZ-BG _ e e;^4 

^00 



K^dkpe pQ^ 



oo 



A=LA,TA 



dqz / dg- _ 2 

-oo Jo (cQ — Ij 

2^l'r^\rr5^2^2 / OT 29 



i2K-^dkpelql ^ -6 + „6 ) '^^■^•^^ 



with as the Riemann function: ^(5) ~ 

1.037. To derive Eg. (|A.13|) . the (/7q-independent 



forms of the scattering matrix are used: IMqaIpJoz = 
ATT-'e-'el^Bx/ [tedusx\e{q,nQx)\''] [Bla = 9g4g2/Q7 

BtA = (g^ + 8g2g4)/Q7] 41 ^j^j^j^ jjfjgj. gjightly from 

those presented in Sec. III. 

Using the material parameters for a GaAs/AlGaAs 
heterojunction, Qqo^"^*^ and (5o(f ^^'^ vary with temper- 



ature as Qoo^"^" « 1.936 x IQ-^TV^s^^ and Q, 
3.085 X lO-^TVris^^ {ris in 10" cm'^), respectively. 
This implies that, of the various electron-phonon scatter- 
ings, the (longitudinal-phonon) piezoelectric interaction 
is dominant at low temperature. 

In the EP regime, Tbg < T ep and uqx ~ uqx + 
1 w T/f^QA, and Eq. (jA.lOp reduces to 

2kF 



PZ-BG 
00 



T 



inkpe'p 



E 



IH^qz^dq, 



xlM, 



QX\ 



qyn<ix 



1 - 

J- 4p2 



(A.14) 



Accordingly, the contribution to Qqq from defor- 
mation and piezoelectric scatterings in the EP regime. 



DP-EP , ^DP-EP • V, 

and Qqq , are given by 



'00 



and 



^DP-EP 
^00 



\I{iqz)?dq,, (A.15) 



^PZ-EP 
Voo 



K^depkp 



1 f°° 
— / \I{2lkpq,)\''dq, 
^sX J~oo 



dq- 



14 

A 

fBx/Q 



, . (A.16) 

[l + qsl{2kpq)]^,/T^ 

Thus, we find that the electron-acoustic-phonon scat- 
tering tends to p ~ Qqq ~ T^ in the BG regime and 

P ~ Qoo ~ T in the EP regime. 

It should be noted that the resistivity correction we 
found here in the analysis of electron-impurity scatter- 
ing, which goes beyond the earlier "linear-in-T" result, 
is consistent with previous transport studies in Refs.lM 
andlrsl. In regard to electron-phonon scattering, our re- 
sult concerning the power-law temperature dependence 
of resistivity due to piezoelectric scattering agrees with 
the previous one: in Refs. [tI and [tI the contribution 
to inverse relaxation time due to piezoelectric scatter- 
ing was found to be proportional to T^. However, our 
deformation-scattering result is different from the T'' law 
obtained previously. This is associated with the fact that 
the deformation scattering is taken to be unscreened in 
present paper, while a screened one was used in the pre- 
vious studies i^^ii^ 



2. Diffusion thermoelectric power in the 
high-electron-density limit 

To obtain the diffusion TEP Sd, one needs to consider 
the second term on the left hand side of Eq. ([15]). This 
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term can be written for n = as 

-jr- 2^ aom7m = -V^r (^1 - — j e ^ 

m 

_V,T— e-^ +ofe- 
2T V 



and for n = 1 we have 

^-^ 3 TT 



(A.17) 



T 



O (e^ — 

(A.18) 



It is clear that, in high-electron-density limit, the n = 
term is exponentially small while the term with n = 1 is 
dominant and is given by 



T 



E 



(A.19) 



Assuming r]n{ep) ^ rjn{ep-q) in the case T, Tbg < ^f, 
Qj^?'^*^ can be written as (n, n' < 1) 



-\imp,ph 
-' nn' 



poo 

/ de 




'o 


de 



V3{£f ~ e) 



71 + ^?-' 



A'"^P'P'Xe).(A.20) 



Using the low-temperature expansion of the Fermi 
function^^ 

6 30 

(A.21) 

and performing the energy integration, the leading terms 
of Q"^'^^ take the forms 



^imp,ph ^imp.ph 

ho — ^01 



3 dep 



,(A.22) 



and 



imp,pli ^imp,pli ^imp,pli 



\eF). (A.23) 



Substituting these Qnn' terms into Eq. ()15p , the solution 
Cq*^^ can be written as 

Cl,"^ ~ ^-T^T-^ [A'-P(£j.) + AP'^(e;^)] , (A.24) 



and the diffusion TEP takes the form 



3e 

fT_d_ 
3e dep 



dep A™P(eF)lT=o 
T{ep)T + K^\ep) 



A""p(£f)It=o 



(A.25) 



From Eq. (jA.25l) we see that the first term on the right 
hand side agrees with the well-known Mott formulai^ 
Considering the fact that Qoo relates to Us approxi- 
mately as Qoo ~ "T-* {'^^^) /{Qs^s) and ignoring the 
energy-dependence of (t^^), Sd ~ ■k'^T /{ieep) can be 
obtained. However, if one assumes ^ e^, we ob- 

tain 5^ ~ TT^T I (2te£w){l in agreement with the re- 

sults of Refsl45l andl46l. 

In Eq. (|A.25p , the second term on the right hand side 
is a low-temperature correction to the leading term and 
it is proportional to T^. Obviously, in the BG regime, 
this correction comes only from the temperature depen- 
dence of the screening of electron-impurity scattering, 
since the phonon contribution is proportional to in 
BG regime [A^'^{£p) ~ T^] and it can be ignored. How- 
ever, in the equipartition regime, the electron-phonon 
scattering results in AP^(ei;') being linear in T. Hence, 
both the electron-impurity and electron-phonon scatter- 
ings lead to a deviation of Sd vs T from the linear rule 
when Tbg < T < e^- 



3. Phonon-drag thermoelectric power in the 
high-electron-density limit 

To investigate the phonon-drag effect in thermoelec- 
tric power, one needs to study the driving term Z?„ in 
Eq. (fT5|) . Performing substitution, q q for + term 
and p — > p + q for — term, Eq. (ITCl) can be rewritten as 

Ai = 27r ^ |MQAp(5(£p+q - £p - ^Q\)Tp\UQx 



T 

X [Xn{p 



nQ\{l + nqx) [/o(ep) - /o(£p+q)] 

q)-Xn{p)]. (A.26) 



Considering only the driving term £>„ with n = 0,1, 
the solution of Eq. P3|) . Cq^\ can be written as 



Oq — 



T QiiQoo 
V,T f 1 



T 



QoiQio 
Qoi 



Qoo " QooQii 



(A.27) 



Using Eqs. (|A.22I) and (|A.23|) . the phonon-drag Sg is de- 
termined by 



Sn 



1 



Do + ^^TD^^\nK{£p) 

3 OEp 



(A.28) 



Recognizing that Xn{p + 9) ~ Xn{p) ~ 
Attq cos ipqrin{ep)/p'^ in the case T <^ ep, Sg finally 
takes the form 



27r^ dlnAjep) ^ ~ 



de 
QT 



Q,P,X 



nQx{l + UQx) [/o(ep) - /o(ejj(Afl9) 
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Note that this expression for Sg reduces to the one widely 

used in hterature^^i^iS if is replaced by l/ep- 

To further analyze the power faw of Sg vs T in the high- 
rig limit, one has to study the temperature dependence 
of TpA- At sufficiently low temperature, it is reasonable 
to assume that boundary scattering dominates phonon 
relaxation and the mean free path of phonons. A, is in- 
dependent of T. Under this consideration, Sg in the BG 
regime, S^'~^, can be written as 



S, 



BG 



-OO Jo 



dq 



X \Mqx\ 



2 Tpxq'^Qe^ 



(A.30) 



:.EP 



and, in the EP regime, it takes the form 
: \MqxY 



Sg can be further simplified by substituting explicit 
forms of the deformation and piezoelectric scattering ma- 
trices into it and then performing momentum integration: 



S. 



BG 



15m*^S^rVp,LAC(5) 



From Eqs. (|X30)) and (|X3T|) we see that, when T ^ 0, 
the phonon-drag thermoelectric power tends to zero as T"^ 
in the BG regime and it reaches a saturation value in the 
EP regime. Note that such behavior of Sg vs T in the BG 
regime has already been demonstrated in RefsIH andlTSl. 
while, as far as we know, the temperature-independence 
(A. 31) of Sg in the EP regime, obtained here, is a new prediction. 
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